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Abstract 

This article discusses the possible variation of the solar neutrino flux over time in the 
records of Super-Kamiokande-I and the relation to non-equilibrium statistical mechanics the 
and entropic pathway model. The scaling behavior of the Super-Kamiokande-I time series is 
investigated utilizing Standard Deviation Analysis and Diffusion Entropy Analysis. The data 
set exhibit scaling behavior and may follow Levy-type probability distribution function. 
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1. Introduction 

In recent years, physicists and mathematicians have been discussing the possible signature 
of variation of the solar neutrino flux over time in the records of solar neutrino experiments 
and the relation of such a signature, if discovered, to non-equilibrium statistical mechanics 
as developed by Tsallis [TO] and to the entropic pathway as introduced by Mathai [1]. 
Spatio-temporal analysis of solar activity and possible correlation to the variation of the solar 
neutrino flux, as reported in the literature, have also been debated [11 J. Such discussions 
were central to the UN/ES A/NASA/ J AXA workshop for the preparation of the International 
Heliophysical Year 2007, Abu Dhabi and Al-Ain, United Arab Emirates, 20-23 November 
2005 [9J. This article reports some progress in the discussions made since then. 

2. Time Series: Variance Versus Probability Density Function 

Time series are recording the evolution of a parameter over time and commonly exhibit 
stochastic fluctuations of the recorded parameter if the underlying process is random. Tradi- 
tional methods for the study of scaling exponents of such time series like standard deviation 
analysis and spectral analysis (Fourier, wavelet) are based upon the analysis of the variance 
of the time series [SIE]. 

The statistical variability of the average solar neutrino flux is interpreted as noise (see 
Fig.l). The solar neutrino flux fluctuations are thought to contain no useful information and 
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are consequently smoothed and show a constant average in the long time. According to the 
central limit theorem, the statistics of the fluctuations in such a time series, emanating from 
complex systems such as the gravitationally stabilized solar fusion reactor, are assumed to 
be Gaussian. The fact that they may not be, due to the fact of being generated by an open 
system, remains unexplained. However, suspected non-Gaussian behavior prompted studies 
of possible time variation of solar neutrino records as a problem belonging to non-equilibrium 
statistical mechanics wherein statistical fluctuations can provide useful information about, 
for example, the reaction and diffusion properties of complex phenomena. An example would 
be the fluctuation-dissipation theorem, in which the response of a complex system to a per- 
turbation is determined by the complex system's unperturbed autocorrelation function. The 
variation of the solar neutrino flux would be indicative of the Sun's complex internal dynamics 
as evidenced by changes in the number of detected solar neutrinos. Such studies maintain 
that the variations in the solar neutrino flux detected on Earth are not noise, but contain 
information about the source of variability, particularly the variations in the gravitationally 
stabilized solar fusion reactor, beyond neutrino oscillations [8J. 

Different from standard deviation analysis (SDA), Scafetta's method of diffusion entropy 
analysis (DEA) evaluates, under certain conditions, the scaling exponent of the probability 
density function (pdf) through the Shannon entropy (or Tsallis entropy, or Mathai's path- 
way) of the diffusion process generated by the fluctuations exhibited by the time series [7], 
such as shown in Fig.l. The application of DEA shows that the fluctuations in the solar 
neutrino record may have inverse power-law statistical distributions. Specifically, if x is the 
solar neutrino flux, the distribution of the values of neutrino fluxes is an inverse power law 
P{x) ^ where C is a normalization constant. The parameter a is the inverse power-law 
index. The scaling of the solar neutrino record was tested by randomly changing the order 
of the data points. If the time series were internally correlated, the resulting distribution 
would have changed from the original, but that could not be observed. The invariance of the 
distribution under shuffling indicates that the statistics of the time series in non-Poisson and 
renewal, meaning that with the generation of each new event, the process is renewed [12j . 

3. Diffusion: Gauss Versus Levy Statistics 

The statistics of the solar neutrino record is described by a non-Gaussian distribution. 
The behavior of such limit distributions requires a generalization of the central limit theorem 
to the case in which the second moment of the variate diverges. Such processes were studied 
by Levy. The solar neutrino record is describable by such a Levy distribution. 

The diffusion process, as utilized in Scafetta's DEA, is modeled as a continuous-time 
random walk where a random walker jumps instantaneously from one location in space to 
another following a waiting period on a location whose duration is drawn from a pdf of 
waiting times [5, [l^. Depending on the rules of walking in space and jumping over time, 
pdf's of the diffusion process show Gaussian behavior or can be represented by Levy-type 
flights or Levy- type walks [3]. Such pdf's can be derived from Mathai's entropic pathway. In 
Mathai's pathway model, a specific parameter allows to proceed from probability densities of 
a generalized type-1 beta model to a generalized type-2 beta model to a generalized gamma 
model when the variable is restricted to be positive [4J- See Figures 1 and 2. 
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The Figure 1 (a), (b) and (c) show the behaviour of the generahzed type-1 beta form of 
the pathway model given by the function fi{x) = cix^[l — a(l — a)x ] i-« for the values 
7 = 2, a = 2, 5 = 2 and for a = 0,q: = 0.5 and a = 0.8 respectively, ci is the normalizing 
constant. This is a model with the right tail cut-off when a < 1. 
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(a) - (b) 

The Figure 2 (a) shows the behaviour of the generalized type-2 beta form of the pathway model 
given by the function f2{x) = C2X^[l-\-a{a — l)x^]~^^ for the same values 7 = 2, a = 2, 5 = 2 
and for a = 1.1, a = 1.2 and a — 1.4, respectively. Here C2 is the normalizing constant. 
Superstatistics of Beck and Cohen [I] belongs to this case, (b) shows the behavior of the 
generalized gamma form of the pathway model or the Gaussian model given by the density 
function /3(x) = C3X^e~^^ for the values 7 = 2, a = 2, 5 = 2. As o: ^ 1 the type-1 and type-2 
beta models fi{x) and f2{x) reduce to Gaussian form given by f^{x). Taking into account 
correlations (memory), deviation from Gaussian behavior can be described by application of 
fractional derivatives in space and time. The diffusion process can be extended to a reaction- 
diffusion process if the species are not only mobile in space and time but are also allowed to 
react with each other [4J. In this case, diffusion, i.e. random walk of reactants, is the mecha- 
nism that allows for reaction, in more general terms described by fractional reaction-diffusion 
equations. Closed-form representations of such equations, reflecting non-Gaussian behavior, 
are achieved by utilizing Mittag-Leffler functions and H-functions [2, 4J. 
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Figure 3 shows the horizontal axis with a time period (1997 — 2002) from the beginning of the 
data taking and the vertical axis which is the measured neutrino flux in units of 10^cm~^s~^. 
The correction is included in the measured solar neutrino flux of 5-day-long samples from 
Super-Kamiokande-I [13j. It is commonly assumed that this record of data does neither 
contain a trend nor a periodicity. The average value of the record confirms the expected 
average solar neutrino flux based on a standard solar model. 
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Flux Variation - SDA 
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The Figure 4 shows the results of standard deviation analysis (SDA). In the SDA (blue line) 
one examines the scaling properties of the second moment of the diffusion process generated 
by the time series as shown in Fig. 1. The standard deviation D{t) of the variable x is 

D{t) = (< x^t > - < >2)^ - t^. 

The Hurst exponent H is interpreted as the scaling exponent being evaluated from the gradient 
of the fitted straight line in the log-log plot of D{t) over diffusion time t. In the Figure, the 
Hurst exponent oi H — 0.58 for the red straight line does indicate a deviation from Gauss. 
Also shown in green is the line for H — 0.5. 
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Flux Variation - DEA 

1 > > > > — > — 



DEA 




Time (log) 



The Figure 5 shows the results of diffusion entropy analysis (DAS). The DEA determines the 
scaling exponent delta evaluated through the Shannon entropy S(t) of the diffusion process 
generated by the variation of the time series in Fig. 1. The probability distribution function 
t) is determined by means of the sub-trajectories 

/ 

^n{t) = ^Ci+n, n = 0, 1, . . . . 
1=0 

If the scaling condition 

is valid, the corresponding entropy increases over time as 

r+oo 

S{t) = - / dx p{x,t) ln[p{x,t)]. 

J oo 

Using again the scaling condition for p{x, t) one obtains 

S{t) ^A + S\n{t). 
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With 

r+oo 

A = / dy F{y) ln[F{y)] = constant, 

J oo 

where y = fs- This indicates that in the case of a diffusion process with a scahng proba- 
bihty distribution, its entropy S{t) increases hnearly with ln(t). The scahng exponent 5 is 
being determined from the gradient of the fitted straight hne in the hnear-log plot of S{t) 
over time t. For fractional Brownian motion the scaling exponent 5 coincides with the Hurst 
exponent H. For random noise with finite variance, the probability distribution function 
p{x^t) will converge, according to the central limit theorem, to a Gaussian distribution with 
H = 6 — 0.5 . If the Hurst exponent H is different from the scaling exponent 5, the scaling 
represents anomalous behavior. The diffusion process characterized by Levy flights belongs 
to anomalous diffusion. In this specific case the scaling condition for p{x^ t) is still valid but 
the variance is not finite and the variance scaling exponent can not be defined. A second 
type of anomalous diffusion is due to Levy walk based on a generalization of the central limit 
theorem. In this case of anomalous diffusion, the second moment is finite and the exponents 
H and delta obey the relation 5 = 1/(3 — 2H). In the Figure, the scaling exponent oi 5 = 0.52 
for the red straight line does indicate a deviation from Gauss. Also shown in green is the line 
for 5 = 0.5. 
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